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Abstract 

The phenomenological dissipation of the Bloch equations is reexamined in the context of 
completely positive maps. Such maps occur if the dissipation arises from a reduction of a unitary 
evolution of a system coupled to a reservoir. In such a case the reduced dynamics for the system 
alone will always yield completely positive maps of the density operator. We show that, for 
Markovian Bloch maps, the requirement of complete positivity imposes some Bloch inequalities 
on the phenomenological damping constants. For non-Mar kovian Bloch maps some kind of 
Bloch inequalities involving eigenvalues of the damping basis can be established as well. As an 
illustration of these general properties we use the depolarizing channel with white and colored 
stochastic noise. 

PACS numbers: 42.50.Dv, 03.65.Ud, 42.65.Lm 
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I. INTRODUCTION 



In 1946, Felix Bloch introduced a set of equations describing the dynamics of a nuclear 
induction of a spin that interacts with a magnetic field Qj. The applied magnetic field drives 
the Bloch vector of the magnetic moment, causing it to precess about the field direction. In 
addition to the unitary evolution describing the magnetic moment precession, a nonunitary 
evolution is observed in nuclear magnetic resonance, which results in dissipation of the 
magnetic observables. This dissipation is characterized by two phenomenological decay 
constants. The two lifetimes 7\ and T 2 are the longitudinal and transverse decay constants, 
respectively. 

Fluctuations in the environment, such as inhomogeneities in the magnetic field and in- 
teractions with other moments, lead to dissipation in the system. The dissipation constants 
Ti and T 2 are nonnegative so that exponential decay of the magnetic moment occurs. It is 
well known that this condition is required to preserve the positivity of the density operator 
under dissipation. The phenomenological dissipation of the Bloch vector defines a positive 
map (PM) of the density operator for the spin system. 

About the physical sources of the dissipation constants Bloch wrote: 

The actual value of T\ is very difficult to predict for a given substance ... To 
give a reliable estimate of T 2 ... requires a more detailed investigation of the 
mechanism involved and will not be attempted here. 

In the same paper Bloch made the following statement about the relative values of the 
dissipation constants: 

... serious errors may be committed by assuming 7\ = T 2 . There are, on the 
other hand, also cases where this equality is justified ... 

It took almost 30 years to understand that the dissipation results from a reduction of 
a unitary evolution of a system coupled to a quantum reservoir. In the process of such a 
reduction, the transformation of the density operator of the system has to be a completely 
positive map (CPM). However, if the hypothesis of complete positivity is to be imposed 
then the values of the dissipation constants cannot be arbitrary. In particular, the inequality 
2Ti > T 2 must hold and has been experimentally observed Q|. 
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It is the purpose of this paper to reexamine the well-known Bloch equations in the con- 
text of completely positive maps. We show that the condition of complete positivity for 
Markovian Bloch maps imposes some Bloch inequalities on the phenomenological damping 
constants. For non-Markovian Bloch maps, generalized Bloch inequalities involving eigen- 
values of the damping basis can be established as well. The depolarizing channel with white 
noise is used to illustrate these general properties. The non-Markovian Bloch map is studied 
in the framework of a depolarizing channel with colored noise. 

II. BLOCH EQUATIONS 

Although originally introduced in the context of nuclear magnetic resonance, the Bloch 
equations are well-known in quantum optics, where they describe a two-level atom interact- 
ing with an electromagnetic field (3J. The Bloch equations offer a physical picture of the 
density operator. The dynamics of any two-level quantum system can be expressed in terms 
of a three-dimensional vector b = (u,v,w) T , called the Bloch vector. For such systems, the 
set of all density operators can be geometrically represented by a sphere with unit radius. 
States that are on the surface of this Bloch sphere are pure states or rank one density oper- 
ators p = States that are within the Bloch ball are mixed states, which are written 
as convex combinations of pure states. 

The optical Bloch equations are a set of differential equations, one for each component 
of the Bloch vector, having the form 

1 

u = -—u - Av, 

J- u 

v = -^v + Au + Qw, (1) 
w = -—(w-w eq )-fl,v. 

The unitary part of the evolution is governed by Q, the Rabi frequency of the applied field. 
The field is detuned from the natural resonance of the atom by an amount A. Note that these 
equations differ from the original Bloch equations by the fact that there are two different 
transverse dampings. These two constants T u and T v are the decay rates of the in phase and 
out of phase quadratures of the atomic dipole moment, while T w is the decay rate of the 
atomic inversion into an equilibrium state w eq . The interaction Hamiltonian of the system 
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and reservoir that leads to Eqs. (0) is 

H = h{Qa + fiV) + H[aT(t) + ^T\t)}, (2) 

where a and are the lowering and raising operators for the atomic system. The master 
equation for the entire system (S) and reservoir (R) is given by the von Neumann equation 
(ft = l) 

Psr = -^H.Psr}- (3) 

The master equation for the system alone is obtained by tracing over the environment degrees 
of freedom. The decay constants arise due to an interaction with the reservoir given by the 
variables T(t). This could, for example, be a collection of harmonic oscillators, in which case 
T(t) = Hk9k0k e ~^ u ~ Vk ^ t ■ m contrast to a quantum reservoir, T(t) could describe a classically 
fluctuating environment. 

Typically, the phenomenological decay rates in the Bloch equations appear as 

1 = 1 1 = 1 1 = 1 (4) 

u J-2 J-v J-2 J-w J l 

The damping of the component that is in phase with the driving field is equal to the damp- 
ing of the component that is out of phase with the driving field. The damping is caused 
by the interaction of the system with an external environment that is averaged over. This 
environment could be the vacuum field, which leads to spontaneous emission. Other phe- 
nomena, which result in Eq. (Jl|). are coupling to a thermal field or phase randomization 
due to atomic collisions. The source of noise and dissipation in the system is due to the 
fluctuation of the environment. 

The phenomenological decay rates are not limited to Eq. (@J. Reservoirs that result 
in different dynamics can be engineered. For example, a two-level atom interacting with 
a squeezed vacuum reservoir will experience unequal damping for the in phase and out of 
phase quadratures of the atomic dipole. The corresponding damping rates become 

11111111 

T, rji rji 5 rji rji rri 5 rri rji ' V / 

u ±2 J 3 J-v ±2 J 3 J-w J l 

The presence of the parameter T 3 is the source of the damping asymmetry between the u 
and v components of the Bloch vector. It arises because the vacuum is squeezed, meaning 
that it has fluctuations in one quadrature smaller than allowed by the uncertainty principle 
at the expense of larger fluctuations in the other quadrature. Equation (0) reflects this - 



the decay rate for one component of the Bloch vector is increased, while the decay rate for 
the orthogonal component is correspondingly decreased. The physical parameters leading 
to this dynamics are 



i-^+l+IM,), ±^(iV + i-|M|), (6) 



' 2 A ( N + 1 



T w \ 2 J ' H 2N + 1 

where A is the Einstein coefficient for spontaneous emission, N is the mean photon number 
of the squeezed vacuum reservoir, and M is the amount of squeezing of the reservoir. The 
parameter N is related to the two-time correlation function for the noise operators of the 
reservoir (a/(t)a(t')) = NS(t — t') where a(t) is the field amplitude for a reservoir mode. The 
squeezing complex parameter M arises from the two-time correlation function involving the 
square of the field amplitudes (a(t)a(t')) = M*5(t — t'). The damping asymmetry parameter 

= A\M\ is in this case due entirely to squeezing. These are the relations obeyed by 
squeezed white noise, which lead to squeezing of a vacuum reservoir j^. 

The solution to the Bloch equations determines the state of the system for all time. The 
Bloch vector evolves in time according to a linear map. This linear map may be written in 
the form 

$:bt-+P = Ajb + t, (7) 

where A is a damping matrix and t is a translation. The overall operation consists of 
contractions and translations. Due to the presence of translations, the transformation is 
affine. 

The damping matrix is a 3 x 3 matrix that takes the diagonal form 



A 



( Ai ^ 

A 2 
\0 A 3 j 



(8) 



We will show in Section V, that these eigenvalues can be calculated using the damping basis 
for an appropriate master equation describing the unitary and the dissipative dynamics of 
the spin system. 

Due to the correspondence between the Bloch vector b and the density operator p, the 
linear map is a superoperator that maps density operators into density operators according 
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to 

$ : p i-f p' . (9) 

The density operator can be expanded in the Pauli basis o a = {I, a x , a y , a z } and the com- 
ponents of p transform under the map. This transformation is characterized by a 4 x 4 
matrix representation of $. It has been found that the general form of any stochastic map 
on the set o: 
parameters 

Without loss of generality, this 4x4 matrix may be cast into the form 



complex 2x2 matrices may be represented by a 4 x 4 matrix containing 12 

30 




(10) 

which uniquely determines the map. Hermiticity of the density operator is preserved by 
requiring that T be real. The first row must be {1,0,0,0} to preserve the trace of the 
density operator. 

The matrix representation of the Bloch vector as an expansion in terms of the Pauli 
matrices 

/ w u — iv \ 
B = b.a=\ (11) 

\ u + iv —w I 

illustrates the properties required by the linear map. In the absence of noise, the Bloch 
vector remains on the Bloch sphere so that 

detB = -(u 2 + v 2 + w 2 ) (12) 

has magnitude unity. A general map transforms the matrix B according to 

<&:B^B'. (13) 

To guarantee that the map $ transforms the density operator into another density operator, 
the Bloch vector can be transformed only into a vector contained in the interior of the Bloch 
sphere, or the Bloch ball. This requirement implies 



|detS'| < |detS|, (14) 



so that the qubit density operator 



p= l -(l + b.a) = \{I + B) (15) 



under the map becomes 

$(p):p-$(p) = ± (/ + £')■ (16) 
This is only possible if the Aj in the damping matrix contain contractions. This is achieved 
in general for |Aj| < 1. If Aj = e~*/ Ti then 1/T, > for i = u, v, w is necessary for $ to be a 
positive map; i.e., it always maps positive operators into positive operators. 

The set of all pure states lie on the surface of the Bloch sphere u 2 + v 2 + w 2 = 1. The 
map $ takes this set into a set of states that lie on the surface of an ellipsoid 

Thus, it typically maps pure states into mixed states. All ellipsoids that are on or inside the 
Bloch sphere represent sets of positive operators. However, not all ellipsoids on or inside the 
Bloch sphere correspond to completely positive dynamics. 



III. COMPLETELY POSITIVE MAPS 

n 

The map $ should be a completely positive map |7(. A completely positive map is defined 

by 

$®/ n >0VneZ+, (18) 

where the index n is a positive integer in the set of all positive integers Z + . The definition 
states that if the dynamics $ occurs on the system and external systems are attached to the 
system, which evolve according to the identity superoperator /, and if the the overall state 
is positive then the map $ is completely positive. 

A completely positive map is required to describe reduced dynamics because it implies 
that the reduced dynamics arises from a unitary evolution 

%{p) = Tr r {U(t)(p ® |7o)(7o|)f/ t W}, (19) 

on the larger Hilbert space consisting of the system and the environment. The environment 
degrees of freedom are denoted by V and 70 is some initial state of the environment. Starting 
with a Hamiltonian for the closed system plus environment, and then tracing or averaging 
over the environment degrees of freedom, will always yield completely positive reduced 
dynamics for the system alone. Because of this requirement, there are points inside the 
Bloch sphere that are not accessible. 
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If the reduced dynamics is consistent with Eq. (jl9j) then it also has as a Kraus decom- 
position. This implies the existence of a set of operators Ki, called Kraus operators, such 
that the map can be expressed as 

<$>{p)=Y,K\ P K u (20) 

i 

where the condition 

Y^K t Kj = I (21) 

n 

ensures that unit trace is preserved for all time p|. If an operation has a Kraus decomposi- 
tion, then it is completely positive. The converse is also true. 

To check whether a map $ that takes n x n matrices into n x n matrices is completely 
positive, it is necessary and sufficient to check the positivity on a maximally entangled 
n 2 x n 2 state ^. This is a powerful theorem that provides a test on a finite space, rather 
than relying on the less practical definition, which requires attaching systems in a countably 
infinite space. In addition, it makes no reference to Kraus operators, but rather, guarantees 
their existence. 



IV. BLOCH INEQUALITES 

A general completely positive, trace-preserving map for two-level systems can always be 
written using four or fewer Kraus operators. An important class of maps called unital maps 
(no translations) has the following set of Kraus operators: 

K = ^ v /l + Ai + A 2 + A 3 I 

K x = l y /l + A 1 - A 2 - A 3 o x 

\ i (22) 

K 2 = -v 7 !- Ax + A 2 -A 3 a 2 



ir 3 = -V / l-Ai-A 2 + A3 o- 3 . 



S 



In order that the class of unital maps be completely positive, it must be that the expressions 
under the radical are nonnegative. The damping eigenvalues must obey the four inequalities 



A!+A 2 -A 3 < 1, 



Ai -A2 + A3 < 1, 



(23) 



-Ax + A 2 + A 3 < 1, 



— Ai - A 2 - A 3 < 1, 



to guarantee complete positivity of the map. This is a necessary condition. This condition is 
more restrictive if compared with the condition: |Aj| < 1 required by a positive map leading 
to contraction of the Bloch vector. In Figure 1 we have depicted the completely positive 
maps as points inside of a tetrahedron, forming a subset of all positive maps contained in a 
cube. We shall call the relations (J23J) Bloch inequalities. 



FIG. 1: Geometric representation of PM versus CPM for the Bloch equations. All points of the 
unit cube represent positive maps. Points inside the tetrahedron are completely positive maps, 
given by the conditions (|23jl . 

The Bloch inequalities (|23J) lead to certain restrictions on the damping constants in the 
Bloch equations. For purely exponential character of the A i; these conditions lead to a 
simpler relation involving only lifetimes. In this case each component of the Bloch vector 
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can only decay according to 



1 1 

T ~~ T ' T 



1 1 

Ti — rri rri 

W -*-u J-v 

11 1 



(24) 



The condition of complete positivity leads to a set of Bloch inequalities for the phenomeno- 
logical lifetimes that must be satisfied. This explains the well-known phenomenon in nuclear 
magnetic resonance whereby the inverse transverse relaxation time is always less than or 
equal to twice the inverse longitudinal relaxation time. 

It is instructive to see how the violation of the Bloch inequalities (|2*1|) is reflected in the 
general conditions (j23j). In order to show this, we have depicted in Fig. 2 the four conditions 
()23|) as a function of time for two different selections of the damping parameters. A function 
exceeding the straight line at the value one is not a CPM, because it violates the Bloch 
inequalities. Note that the case describe on the left figure is a PM that is not a CPM for 
all times, while the case on the right figure describes a CPM for all times. Application of 





FIG. 2: Plots of the four conditions Q23|) as a function of time t. The left figure corresponds to 
7tt = 6, = 3, Jh- = 1. The right figure corresponds to J- = 6, Jr = 5, J- = 1. All values 



arc 



in arbitrary units. We see that the figure on the left is not a CMP, because the Bloch inequalities 
J2U) are violated. 



the Bloch inequalities to the case of squeezed noise given by Eqs. (JHJ) leads to the condition 



N+ \ > \M\ 
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As a different example of a Bloch CMP generated by a noise, let us consider the Hamil- 
tonian 

H = x{t)a x + y(t)a y + z(t)a z , (25) 
where x(t), y(t), and z(t) are Gaussian random variables with 

r„ 



(x(t)x(a)) = -fS(t - s) 



(y(t)y(s)) 

(z(t)z(s)} 



2 

r, 



S(t - s) 
6(t - s) 



(26) 



Note that the white noise correlation has the form (a(t)a(s)) = 2D a 5(t — s), where the 
diffusion coefficient is D a = The same holds for the remaining noises. 

After averaging over the random variables, the Bloch vector will transform by exponen- 
tiating the following matrix: 



/ 



V 



o 
o 



o o 
r a + r c o 
o r a + r b 



\ 



J 



(27) 



As a result we obtain three different lifetimes given by the following relations: 



^ = r fc + r c , l = r a + r 



i 

-*- in 



r a + r 6 . 



(28) 



It is easy to verify that for positive diffusion coefficients these lifetimes always satisfy the 
Bloch inequalities (|24|) i.e., are always generating completely positive maps of the Bloch 
vector. An analysis of other cases in which there are three phenomenological decay rates 
can be found in Ref. 



V. MARKOVIAN BLOCH MAP 

The Bloch equations with Gaussian white noise discussed in the previous section are just 
an example of a general class of Markovian completely positive maps. The general quantum 
Markovian master equation for the system density operator can be written in terms of the 
Gorini-Kossakowski-Sudarshan generator 



n density c 



Lp = -i[H,p]+£p, (29) 
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where the Lindblad superoperator may be written as 



C P = \ll C ^ PF}} + [Ftp, F/]} (30) 
».i=i 

for p G -M n (*M n denotes the set of n x n complex matrices), where if = if Tr{H} = 
0, Tr(Fj) = 0,Tt{F^ Fj} = Sij, and (c^) is a complex positive semidefinite matrix. 

The expression in Eq. (|29jl has been shown to generate a dynamical semigroup ^| that 
has the following properties: 

(i) ||$ t p|| 1 = ||p|| 1 VpGy i + (H), t > 0, 

(ii) $ 1 ®/ n >0VneZ + , (31) 
(to) limtio $ t = i, 

(zv) $ 4 $ s = $ t+a , t,s>0. 

Property (i) states that the map is trace-preserving for all positive operators, which form a 
positive cone V-^iTi). || ■ ||i denotes the trace norm in the space of linear operators on the 
Hilbert space Ti. The definition of complete positivity is given by property (ii). We have 
already encountered this property in Eq. ()18j) . The third property is a statement regarding 
continuity. The map $ is continuous from above and approaches the identity superoperator. 
Property (iv) is the essence of the semigroup property. It states that applying the map from 
time to time s and then applying the map from time s to time t is equivalent to applying 
the map from time to time t + s. Obviously, exponential functions have this property. 

In the following we will show that such Markovian maps can be conveniently written in 
terms of a damping basis. The map $ can be written in the general form 

$(.) = ^A J Tr{Lj.}i? i , (32) 

i . . 

where Lj and i?, are left and right eigenoperators, which form a damping basis [lj|. Lj is 
the operator dual to Ri that satisfies the duality relation 

TriLiRj} = S^. (33) 

This is a complete, orthogonal basis with which to expand the density operator at any time. 

This basis is obtained by finding the eigenoperators of the eigenvalue equation. The right 
eigenoperators satisfy 

CRi = RiX, (34) 
12 



while the left eigenoperators satisfy the dual eigenvalue equation 

UC = XL t . (35) 

Both Li and Ri have the same eigenvalues. 

If we use the eigenoperators Ri of Eq. (|3~3j) with corresponding eigenvalues \, then once 
the initial state is known 

p(0) = ^Tr{L l p(0)} J R J , (36) 

i 

the state of the system at any later time can be found through 

p(t) = e ct p(0) = Tr{L lP (0)}A^ (37) 

i 

= ^Tr{ J R l p(0)}A l L i , 

i 

where A, are functions of both time and the eigenvalues Aj. An arbitrary positive map must 
have the contractions |A$| < 1. 

This method is a simple way of finding the density operator for a given £ for all times. 
The solution of the left and right eigenvalue equations yields a set of eigenvalues and eigen- 
solutions: {Xi, Rj}. Once the damping basis is obtained, it can be used to expand the 
density operator. 

As an example let us investigate a stochastic model of a depolarizing channel with white 
noise. As a simple example, we will consider Bloch Eqs.Q that have the following phe- 
nomenological decay constants: 

1 = 0, ^r=D, 1 = 0, w eq = 0. (38) 

J-l J-2 -1-3 

This describes the process of phase randomization of the atomic dipole caused by atomic 
collisions. This leads to equal damping for the u and v components of the Bloch vector with 
contractions that depend on the parameter D due to collisions. Setting both Q and A equal 
to zero, the Lindbladian for this system is 

C* = -D[a 3 ,[a 3 ,*]]. (39) 

C is the generator of the dissipative dynamics 

$(p) = e tc p(0) (40) 
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and always generates a completely positive dynamical semigroup. The decay constant is 
a scalar parameter that arises after tracing or averaging over the environment degrees of 
freedom. The environment could be a quantum reservoir or a classically fluctuating external 
system. The von Neumann equation 

p = -iz(t)[a 3 ,p\ (41) 

for system and environment, containing a stochastic variable z(t), can be averaged over the 
noise exactly to obtain Eq. ()39j) . The noise leads to bit flip errors in the system. 
The damping basis for the bit flip Lindblad equation is given in Pauli basis by 

R a = -^=a a . (42) 

The left eigenoperators are identical to the right eigenoperators because the Pauli operators 
are self-dual. The damping eigenvalues are given by 

A = 0, Ai = -4D, A 2 = -4D, A 3 = 0. (43) 

The time dependent functions Aj become 

A = l, Ai = e- 4Dt , A 2 = e- 4Dt , A 3 = I. (44) 

The reduced dynamics describe dissipation of the system due to the coupling of it with the 
environment. The dissipation is in the form of pure damping. In the Bloch equations, this 
leads to phenomenological decay constants 

111 11 

- = - = - = 4D, ^ = ^ = (45) 

for each component of the Bloch vector. The inversion, given by the w component undergoes 
no damping. The two orthogonal components u and v undergo equal damping. 

The bit flip error equation results in a completely positive, trace-preserving map. Thus, 
there exists a set of Kraus operators, which can be used to write the map as a decomposition 
in terms of these operators. Two Kraus operators can be used to write the decomposition 
as 

$(p) = K oP Kl + K 3 pKi (46) 
with Kraus operators explicitly given by 



' 1 _|_ P -iDt 1 _ p-ADt 

Ko = \l /, K, = ^— a 3 . (47) 
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It is easy to show that these operators are normalized 

K\K + KlK 3 = I, (48) 
thus producing a trace-preserving completely positive map. 

VI. NON-MARKOVIAN BLOCH MAP 

The damping basis method can be applied to non-Markovian Bloch maps. Such maps 
occur for example in a depolarizing channel with colored noise. A bit flip error equation that 
does not rely on the white noise assumption can be derived. This is achieved by assuming 
the environment has correlations that are not of the form of a delta-function, resulting in 
what is called colored noise. The master equation is no longer of the Lindblad type; instead, 
it contains a memory kernel operator leading to a master equation of the general form 

p = KCp, (49) 

where K is an integral operator that depends on time of the form Kef) = L k(t — t')cj)(t')dt' . 
The kernel function k(t — t') is a well-behaved, continuous function that determines the type 
of memory in the physical problem. The solution to the master equation can be found by 
taking the Laplace transform 

sp(s)-p(0)=K(s)£p(s), (50) 

determining the poles, and inverting the equation in the standard way. 

To illustrate this class of master equations, the exponential kernel function 

k(s,t) = e- t ¥ [ (51) 

is used. Consider the Hamiltonian H = z(t)as with a random telegraph signal (RTS) random 
variable z(t) = a(— l) n W. The random variable n(t) has a Poisson distribution with a mean 
equal to t/2r, while a is an independent coin-flip random variable |l6[. The equation of 
motion for the density operator is given by the commutator 

p=-i[H,p] = -iz(t)[a 3 ,p]. (52) 
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Taking the ensemble average over the random variables z(t) leads to an equation of motion 



for the average density operator 



P 



-a e 
'o 



[^3, [<T 3 ,p(s)]]ds. 



(53) 



The brackets denoting the ensemble average are omitted and it should be understood that 
we are considering the average density operator. The master equation for the bit flip error 
with an exponential memory kernel is exact. 

The damping basis diagonalizes the equation. This leads to a single equation, for the two 
nontrivial components, that has the following form: 



d of !*-»] 

_c(t) = -4a 2 y dse ~c(s)ds. 

Taking a derivative gives a second order differential equation 

1 d 

^(t) + --c(t)+Aa 2 c(t)=0. 

The function c(t) is the solution to a damped harmonic oscillator equation of motion 

sin(0£) 



(54) 



(55) 



c(t) = cos(fit) + 



2Qt 



(56) 



where f2 = y 4a 2 — and c(0) = 1. 

The image of this random telegraph signal (RTS) map is similar to that of a depolarizing 
channel with white noise. The transformation given by Eq. (|10|) for this noise is 

( \ o o o\ 

c(t) 
c(t) 
\0 1/ 

which has the same form as the depolarizing channel. Notice that there are no translations 
of the Bloch vector and 



T 



(57) 



A(t) = Ai = A 2 = c(t), A 2 = c(t), and A 3 = 1. 



(58) 



The depolarizing channel with white noise has simple contractions so that c(t) is an ex- 
ponential function. This is a property of white noise. The RTS channel has colored noise 
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leading to a nonexponential function for A(£). Rather, it contains oscillating terms with an 
exponential envelope. A power series expansion gives 



A(u) = l-2a 2 t 2 + 0(t) 3 , (59) 

which shows that the linear term in t is missing. Thus, the standard white noise diffusion 
term vanishes. This is in contrast to the Markovian case, where the functions are purely 
exponential functions in time with parameters defining the characteristic lifetimes. This is 
a general property of the memory kernel and a fundamental difference between white noise 
and colored noise. The white noise limit can be recovered from (|58|) with a singular limit of 
r — > and 7 = 4a 2 = const. In such a white noise limit A wn = e -7 *. 

In terms of the dimensionless time v = t/2r we can write A(u) = e~"(cos(//i/) + sm ^ ), 
with ji = a/ (4ar) 2 — 1. The function A(v) has two regimes - pure damping and damped 
oscillations. The fluctuation parameter, given by the product ar, determines the behavior of 
the solution. When < ar < 1/4 the solution is described by damping. The frequency /i is 
imaginary with magnitude less than unity. When ar = 1/4 the function A(z/) = e~ u {l — v) 
is unity at the initial time and approaches zero as time approaches infinity. In addition 
to pure damping, damped harmonic oscillations in the interval [—1, +1] exist in the regime 
ar > 1/4. 

The Kraus operators for the RTS channel are similar to those for the depolarizing channel. 
The non-Markovian Bloch map can be defined as 

$(p) = K\pK x + K\pK 2 

1 i (60) 

= -[l + A(t)]p+-[l-A(t)W<T 3 . 

The Bloch sphere evolves into an ellipsoid according to the equation 

\ 2 / \ 2 

U \ IV 



M)J ' VAM,' +w L (61) 

Unlike the exponential damped solution for white noise, the function can take on negative 
values and is bounded between ±1. 



VII. NOISE AND SEPARABILITY 



The maximally entangled Bell state can be used to test if an arbitrary map for qubits 
is completely positive. The initial entangled state for two qubits is given by a positive 
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semidefinite density operator. If only one qubit of the joint state is subjected to noise, the 
joint state must also be described by a positive semidefinite density operator for all time. 
For n-level systems, positivity of the maximally entangled n 2 -level state implies that the 
map for the n-level part of the state is completely positive. 

The maximally entangled two qubit state can also be used to address the issue of sep- 
arability. The partial transposition map is an example of a map that is positive but not 
completely positive. The partial transposition map applied to an entangled n 2 state involves 
performing the transpose operation on one half of the state, while performing the identity 
operation on the other part of the state. Note that the partial transposition map is not 
a continuous map connected to the identity superoperator. Although the initial entangled 
state is described by a positive semidefinite density operator, the resulting transposed state 
can be negative. Using the Peres criterion of the positivity of the partial transpose, one 
can determine when the state becomes separable 18]. A necessary and sufficient condition 
for the output state to be nonseparable is that the partial transpose map be negative [lflj . 
To check the positivity of the partial transpose it suffices to examine the eigenvalues of the 
operator given by 

Hpab) = U®Tb[$(pab)) (62) 

where T# denotes the transpose of the state of Bob's qubit. The dynamical map $ is first 
applied to the maximally entangled state in the following way: <&(pab) = I a ® §b[pab}- 
The transposition map acts on one half of the entangled state and the output state is given 
by ®{p AB ) in Eq. (JE2J). 

The four eigenvalues of the partial transpose matrix for the bit flip error are 

A(t) _ A(t) 



ei = e 2 = -, e 3 = — , e 4 



(63) 



2' 2 ' 2 

^From this relation we conclude that such a state is separable if and only if A(t) = 0. 
This shows that in the limit of Markovian dynamics, the separability is reached only in the 
asymptotic limit of t — > oo. The situation is remarkably different for a non-Markovian map. 
In Figure 3 we have depicted the separability condition as a function of the color of the 
noise, characterized by the parameter r. 
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FIG. 3: Plots of non-Markovian A(t) and r = 4, (solid line), r = 0.3 (dot-line), white noise case 
r = (dash-line). All plots in arbitrary units with 7 = 1. The state becomes separable for 
A(t) = 0. 

VIII. CONCLUSION 

The Bloch equations are well-known to physicists working in the field of nuclear resonance 
or quantum optics. These equations have been widely used in quantum information theory to 
describe the dissipation of various quantum channels for qubits. The original Bloch equations 
have been derived on a purely phenomenological ground, and phenomenological damping 
constants have been introduced. We have shown that the physical consequences of a quantum 
mechanical description of dissipation leads to complete positivity of the Bloch maps. The 
condition that these maps are CPM, implies a set of Bloch equations for the generalized 
eigenvalues of the damping matrix. In the case of pure Markovian dissipations these Bloch 
equations can be reduced to simple inequalities for the various lifetimes characterizing the 
qubit. We have shown that this approach can be extended to non-Markovian dynamics. 
In this case the Bloch inequalities involve time- dependent eigenvalues of the generalized 
damping matrices. We have illustrated our points with examples of Markovian and non- 
Markovian depolarizing channels for qubits. We have shown a fundamental difference in the 
separability properties of correlated qubits in such channels. 
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